clc
clear
L = 10;
N = 1024;
N2 = N / 2;
multi =  2 * pi * L; 
omega = @(j,t) exp(-1i * abs( multi * j)*t);

om = (1:1024);

x =  createPsiK(2,1,0,L,1024).';


for t = 1:0.0001:1000
    next = x.*(omega(om,t));
    plot(abs(ifft(next)))
    drawnow
end